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Abstract 

The association of long gamma-ray bursts (GRBs) with Type Ib/c supernovae im- 
phes that they explode into the winds of their Wolf-Rayet progenitor stars. Although 
the evolution of some GRB afterglows is consistent with expansion into a free wind, 
there is also good evidence for expansion into a constant density medium. The evi- 
dence includes the evolution of X-ray afterglows (when X-rays are below the cooling 
frequency), the evolution of the pre-jet break optical and X-ray afterglow, and the 
sharp turn-on observed for some afterglows. Recent observations of short bursts, 
which are expected to be interacting with a constant density medium, provide a 
check on the standard afterglow model. Although radio observations do not support 
the constant density model for long bursts in some cases, the evidence for constant 
density interaction is strong. The most plausible way to produce such a medium 
around a massive star is to shock the progenitor wind. This requires a smaller termi- 
nation shock than would be expected, possibly due to a high pressure surroundings, 
a high progenitor velocity, or the particular evolution leading to a GRB. However, 
the need for the termination shock near the deceleration radius cannot be plausibly 
accomodated and may indicate that some long bursts have compact binary progen- 
itors and explode directly into the interstellar medium. 

Key words: Gamma-ray bursts, mass loss, supernovae 
PACS: 97.10.Me, 97.60.-s, 98.70.Rz 



1 Introduction 



Observations of the long duration 
gamma-ray bursts (GRBs), which 
are the primary topic here, suggest 
that they are associated with the 
deaths of massive stars. One line 
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of evidence for this has been their 
association with s u pernov ae (SNe). 
Wooslev fc BloomI (120061 ) list 11 



cases that are good candidates for 
SNe associated with bursts. How- 
ever, 2 recent long bursts have shown 
no e vidence for a SN to faint lev- 
els ( iFynbo et al.l . 120061 ). Another 



line of evidence is the apparent as- 
sociation of the sites of long bursts 
with regions of active star formation. 
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Fruchter et al.l (120061 ) found that the 
positions of GRBs on galaxies are 
more concentrated to the brightest 
pixels than are core collapse super- 
novae, suggesting that the GRBs are 
associated with more massive stars 
than are most core collapse super- 
novae. 

These properties imply that the pro- 
genitors of most explosions are Wolf- 
Rayet (WR) stars. Expectations for 
the surroundings of WR stars are 
discussed in Section [2J Normal SNe 
Ib/c (without a GRB connection) 
are also expected to interact with the 
surroundings of WR stars. The prop- 
erties of the interaction, as observed 
at radio and X-ray wavelengths, can 
provide insight into the GRB case 
(Section[3]). The SNe Ib/c that are as- 
sociated with nearby, low luminosity 
GRBs distinguish themselves from 
the normal SNe Ib/c and are dis- 
cussed in Section |H Afterglow emis- 
sion gives us a prime method of de- 
termining the surrounding medium 
and has been used since the discov- 
ery of afterglows to infer the density. 
This issue is examined in the light of 
recent observations in Section O A 
frequent deduction from the analysis 
of afterglows is that the surroundings 
have a constant density. In this pa- 
per, GRB is used to refer to the long 
duration GRBs. The short bursts 
are discussed in relation to the long 
bursts in Section [61 Possible ways 
of producing a constant density sur- 
roundings are discussed in Section 
[71 Absorption lines in optical spec- 
tra provide another possible window 
on the immediate surrounding of 
GRBs and are treated in Section [HI 
The various issues related to the sur- 



roundings of GRBs are summarized 
in Section [91 



2 The Surroundings of Wolf- 
Rayet Stars 



In this section, I consider the sur- 
rounding medium created by the free 
wind from a WR star, as it is likely to 
provide the immediate environment 
for the GRB. At some point the wind 
is expected to transition to a region 
that results from interaction with 
the surroundings; this possibility will 
be considered in Section [71 Typical 
parameters for a WR star wind are 
a mass loss rate M = 10~^ Mq yr~^ 
and a wind velocity = 10^ km 
s~^. The wind density = Ar'"^, 
where A = M/AnVw, is the crit- 
ical parameter for a high veloc- 
ity interaction and I characterize 
it by A = A/5 X 10^^ g cm~^ = 

(M/10-5 M© yr-^)n^ikmslV2^- 
For G alactic stars, iNugis fc Lamers 
( I2OOOI ) listed mass loss parameters 
for 64 WR stars, yielding an A* range 
of 0.07 to 7.4. The lowest density 
winds are produced by WO stars, 
because of their high values of Vy^, 
up to 5500 km s~^. The s e resu lts are 
based on iNugis et al.l (Il998l ). who 
determined dumping-correcting ra- 
dio mass loss rates, noting that the 
effect of clumping is small at the 
stellar surface, grows to a maximum 
at ~ 5 — lOi?*, and again becomes 
small in the outer wind because of 
the expansion of clumps at the local 
sound speed. 



A difference of the progenitors of 
GRBs with Galactic WR stars is 
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that the GRB progeni tors probably 



have lower metallicity. iModjaz et al. 



koOl. ) found that the nearby GRB /SNe 
are in regions that are systematically 
more metal poor than the regions 
containing core collapse supernovae. 
Metallicities were determined from 
emission line regions close to the ex- 
plosions. The metallicities could be a 
factor of ~ 6 smaller than solar. For 
more distant GRBs, it is not clear 
whether the GRBs have lower metal- 
licities compared to other galaxies at 
a similar redshift z, but metallicities 
Z ~ 0.05 — 0.5 Zq are indicated. In 
recent years, it has been found that 
heavy elements around the Fe peak 
play a role in driving the winds from 
WR stars, so that their mass loss 
rates are Z dependent. In this metal- 
licity range, mass loss rate s from WC 
stars vary as ~ Z~^'^ (jCrowtherl . 
20061), suggesting that values of M 
for GRB progenitor stars are lower 
than the rates for Galactic WR stars 
by a factor of 2 — 3, and the values 
of are lower by a similar factor. 

Another issue is the possible asym- 
metry of the stellar wind. Polar- 
ization studies of Galactic WR 
stars have generally shown an un- 
detectable amount of polarization, 
although ~ 20% show a polarization 
>0.3% that can be interpreted as 
a de nsity contrast of a factor of a 



2-3 (IHarries et al.l . Il998l ). However 



asymmetry could be a significant fac- 
tor for the small percentage of WR 
stars that become GRBs. A plau- 
sible distinguishing feature of the 
GRB WR stars is rapid rotation, so 
that the central core is rapidly rotat- 
ing. For radiation driven winds, the 
higher gravitational acceleration on 



the polar axis can lead to a higher 
radiative flux and mass loss rate on 
this axis, although a lower tempera- 
ture and higher opacity on the equa- 
tor favo r s equ atorial mass loss (e.g., 
Maederl . |2002| ): the density contrast 
from pole to equator can be a factor 



of a few. iMevnet fc Maedeii (120071 ) 



have suggested that higher mass loss 
along the polar axis is needed so that 
the WR progenitor does not lose too 
much angular momentum through 
its wind. 

Overall, the additional effects to con- 
sider for GRB progenitors compared 
to Galactic WR stars do not have 
a substantial effect on the expected 
wind density when compared to the 
large possible range of densities. 



3 Circumstellar Interaction of 
Normal Type Ib/c Super- 
novae 



The interaction of normal SNe Ib/c 
with their surroundings provides an 
interesting case of comparison for the 
GRB case because the driving force 
is generally understood for super- 
novae and there is synchrotron emis- 
sion resulting from the interaction as 
in the the GRB case. In a normal su- 
pernova, the supernova shock wave 
accelerates through the steep density 
profile at the outer edge of the star; 
the acceleration stops when radia- 
tion can stream freely from the star. 
The radiation accelerates the outer 
gas and the radiation dominated 
shock front disappears. The shock 
wave re-forms as a viscous shock in 
the surrounding stellar wind. There 



3 



is a shocked region bounded by a 
reverse shock on the inside and a 
forward shock on the outside. The 
reverse shock is in the steep outer 
power law portion of the supernova 
density profile. The shock fronts are 
plausible sites of particle accelera- 
tion; however, the material entering 
the reverse shock front has a very 
low magnetic field because of the su- 
pernova expansion, so there is some 
question of the efficiency of particle 
acceleration at that site. 

The deceleration of the ejecta by 
the surrounding medium gives rise 
to a Rayleigh- Taylor instability and 
a turbulent region in the shocked 
layer. The magnetic field can be 
built up in this region, although 
it is not clear that the efficiency 
is high. In nu r nerica l simulations, 
Jun fc NormanI (119961 ) found that 
the field is strongest on the smallest 
scales. The energy density in the field 
was limited to ~ 0.3% of the turbu- 
lent energy density, but the result 
was limited by the numerical reso- 
lution. Another possible source of 
magnetic amplification is related to 
instabilities in the coUisionless shock 



waves (e.g., [Bell, I2OOJ). 



The basic hydrodynamic model of 
steep power law ejecta driving an 
interaction shell into a surrounding 
stellar wind can reproduce the ob- 
served radio emission from SNe Ib/c 
if some fraction of the postshock 
energy density goes into relativistic 
particles and magnetic fields, and 
synchrotron self-absorptio n (SSA) is 
impo r tant at early times (IChevalierl . 



19981 : IChevalier fc Franssonl . [20061). 
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Fig. 1. Peak radio luminosity and cor- 
responding age for well-observed core 
collapse super novae. The dashed lines 
give curves of constant expansion veloc- 
ity, assuming synchrotron self-absorp- 
tion at e arly times (upda ted version of 
Fig. 4 in lChevalieil (|l998l ^'). 



In this situation, the peak flux of 



the radio emission gives information 
on the radius, and thus the velocity, 
of the radio emitting region. Figure 
1 shows observed peak fluxes and 
ages of radio supernovae. The dashed 
lines give the velocities of the radio 
emitting regions if the peak is due to 
SSA. If another absorption process, 
such as free-free absorption, is dom- 
inant, the velocity inferred from the 
turn-on is lower than the actual ve- 
locity. In Fig. 1, it can be seen that 
the SNe II have lower inferred veloc- 
ities, which is both because of the 
relatively low velocities in SNe II and 
the importance of free-free absorp- 
tion. There is a large range in peak 
luminosity for SNe II, which can be 
primarily attributed to a range in cir- 
cumstellar density. Multiwavelength 
observations conflrm the large range 
in density. 
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The normal SNe Ib/c show systemat- 
ically higher velocities than the SNe 
II (Fig. 1), which can be attributed 
to some combination of higher peak 
velocities at the time of shock break- 
out, higher mean ejecta velocities 
because of lower ejecta mass, and 
lower deceleration because of lower 
circumstellar density. There is again 
a large range in peak luminosity; 
this may be due to a range in cir- 
cumstellar density, but in this case 
there is no independent evidence 
for a range in density. If the super- 
novae have similar efficiencies for the 
production of synchrotron radiation 
and a range of circumstellar density 
that is comparable to that around 
Galactic WR stars, the observed lu- 
minosity range can approximately 
be reproduced if ~ ~ 0.1 



( 1 Chevalier fc Franssonl . |2006| ) . The 
required magnetic field is high and 
cannot be produced by compres- 
sion of the wind magnetic field (this 
would require a wind energy fiux 
that was completely dominated by 
the magnetic field). 



sion fades because of the low super- 
nova luminosity, and synchrotron 
emission is the most plausible non- 
thermal mechanism. However, an ex- 
trapolation of the radio synchrotron 
emission falls below the observed 
X-ray emission, especially when 
synchrotron cooling of the radiat- 
ing electrons is taken i nto a ccount. 



Chevalier &: Franssonl (120061 ) sug- 



gested a model of particle accelera- 
tion in a cosmic ray dominated shock 
front so that the particle spectrum 
fiattens to high energy. At low en- 
ergy, the particle spectrum is rela- 
tively steep, with energy index p ^ 3, 
in accord with radio observations of 
SNe Ib/c. At high energies, the spec- 
trum becomes fiat. This spectrum 
results in fairly fiat evolution of the 
X-ray emission, while the radio emis- 
sion decreases. 



4 Low luminosity, nearby GRB- 
SNe 



There are a number of X-ray obser- 
vations of SNe Ib/c, but the data 
are much less extensive than at radio 



wavelengths (j Chevalier &: Fransson 



20061 . and references therein). The 
observed luminosities are higher 
than expected from thermal emis- 
sion from interaction with a normal 
WR star wind, so that a nonther- 
mal mechanism is indicated. Near 
maximum light, inverse Compton 
scattering of photospheric photons 
with relativistic electrons is a possi- 



bility ( 


Biornsson & Fransson. 


2004; 


Chevalier & Fransson. 


2006 


. At 



Figure 1 shows that the 3 low lumi- 
nosity GRB/SN events have higher 
velocities of the radio emitting re- 
gions than the normal SNe Ib/c. 
Their positions in the figure suggest 
semi-relativistic velocities. The high 
luminosity nearby event ORB 030329 
had a 5 GHz luminosi ty of 5 x 10^ ° 



erg s ^ Hz ^ on day 10 (iBerger et al. 



20031 ). indicating highly relativistic 
motion in this case. For the lower 
velocity cases, application of the 
synchrotron theory described by 



later times, inverse Compton emis- 



Chevalier fc Franssonl (120061 ) yields 
mass loss densities A*(eB/0.1) of 
0.1 (SN 1998bw), 1.6 (SN 20031w), 
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and 0.02 (SN 2006aj). The theory 
is nonrelativistic, but should still 
yield approximate results consider- 
ing that these objects indicate only 
semi-relativistic motion. The density 
inferred for SN 2006aj is low be- 
cause of the early tu rn-on (see also 



Waxman et all . 120071 ) 



An important issue for the 3 low lu- 
minosity events is whether some of 
the observed phenomena can be ex- 
plained by the supernova, or whether 
a central engine is needed, as in the 
case of normal GRBs. A supernova 
explanation means that emission as- 
sociated with the interaction of the 
fast, outer supernova ejecta can ex- 
plain the observations. Expectations 
for the supernova case depend on the 
ejecta mass and energy for the explo- 
sions. Optical observations of the 3 
supernovae were extensive and there 
are results on the supernova parame- 
ters: 10 Mq and 50 x 10^^ ergs for SN 
1998bw, 13 Mq and 60 x 10^^ ergs for 
SN 20031w, and 2 and 2 x 10 ^^ 



ergs f or SN 2006aj flMazzah et al. 



2006al lbl). The supernova properties 
of SN 2006aj were closer to the nor- 
mal Ic SN 2002ap than to the bright 
SN 1998bw, and its inferred mass 
and energy would be incapable of 
producing the high velocity inferred 
from the radio emission. The impli- 
cation is that the radio emission is 
related to a central engine. The radio 
observations of SN 2002ap suggest a 
low wind density in this case (Fig. 
1), which appears to also apply to 
SN 2006aj. For SN 1998bw and SN 
20031w, the large supernova energy 
allows the possibility of a super- 
nova orig i n for the radio emission. 
Tan et al.l (l200ll ) have discussed such 



a model for SN 1998bw. 

SN 2006aj showed a thermal X-ray 
component over the first few 1000 
sec that has been interprete d by 



rpr 

Campana et al.l (l2006l ) and lWaxman et al, 
(120071 ) as shock breakout emis- 
sion. The temperature was constant 
at ~ 0.17 keV during this time. 
The radiated energy in the ther- 
mal component w as ~ 2 x 10**^ 



(ICampana et all . l2006l : [ij, 120071 ). 



erg s 



This is several orders of magnitude 
larger than would be expected from 
shock breakout from a WR star, as- 
suming no effect of the WR star wind 
flMatzner k McKeel . Il999f ): in addi- 
tion, the duration of the shock break- 
out emission in this case would be de- 
termined by light travel time effects, 
yielding a timescale ~ 10 s, much 



l ess t han observed. ICampana et al. 



(120061 ) addressed this issue by con- 
sidering the progenitor star to be 
surrounded by a dense WR star 
wind, with A^, ^ 20. This wind den- 
sity conflicts with that deduced from 
the radio emission, but the radio 
emission is at later times and it is 
possible that there was a phase of 
dense mass loss just before the ex- 
plosion. With the dense wind, the 
photosphere is formed at r 5x10^^ 
cm. The corresponding light travel 
time, 200 sec, is still less than the 
duration of the thermal compo- 



nent, so ICampana et al.l (120061 ) and 
Waxman et al.l (120071 ) appeal to an 
asymmetric progenitor structure to 
lengthen the timescale. Another is- 
sue is the total energy emitted in the 
thermal component, considering the 
fairly low energy explosion estimated 



for SN 2006aj mentioned above. iLi 



(120071 ) considered supernova shock 
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breakout in a wind and found that 
the energetics presen t a problem for 
SN 20 06ai. However. IWaxman et al. 
( I2OO7I ) attributed the emission to 
the breakout of a mildly relativistic 
shell; it is possible that such shell 
ejection could be generated by a cen- 
tral engine. It thus appears that the 
early X-ray emission from SN 2006aj 
cannot be accounted for by the su- 
pernova and activity of the central 
engine is needed. Whether the ther- 
mal emission can be explained by 
breakout emission remains uncer- 
tain. One issue is whether the con- 
stant temperature emission can be 
produced if the progenitor is highly 
asymmetric. More detailed modeling 
of the emission is needed. 



5 GRB Afterglows 



As discussed in Section [2l the imme- 
diate surroundings of a long GRB 
is expected to be the wind from the 
progenitor star and one would ex- 
pect the afterglow to reflect inter- 
action with such a surroundings. In 
the time before a jet break occurs, 
there are clear differences between 
evolution in a wind medium and in 
a constant density (often referred 
to as ISM for i r iterst ellar medium) 
flChevaher fcTi . l2000h . The cooling 
frequency, where the synchrotron 
cooling time equals the age, increases 
as t^/^ in the wind case, but decreases 
as in the ISM case. This fre- 

quency typically occurs between op- 
tical and X-ray wavelengths, giving 
the expectation that the flux should 
drop more rapidly with time at opti- 



cal wavelenths than in X-rays for the 
wind case. The opposite is true for 
the ISM case. The peak flux, F^m, at 
the typical frequency Um, is lower at 
lower frequencies oc u^^^ in the wind 
case, but is constant in the ISM case. 
This effect can be best observed at 
radio wavelengths because of the 
large range of wavelengths that they 
provide. Finally, the synchrotron 
self-absorption, Ua, drops as t~^/^ in 
the wind case, but is constant in the 
ISM case. Again, radio observations 
are generally needed. 

The application of these differences 
to observed light curves is compli- 
cated by jet effects, which were gen- 
erally found to occur at early times 
(< 3 days) for bursts found dur- 
ing the BeppoSAX era. Models of 
the deceleration of jets have shown 
some features that are not presen t 
in the simple models ( iGranotl . 120071 ). 
so there is uncertainty in the inter- 
pretation. Overall, detailed models 
of afterglows observed during the 
BeppoSAX era generally prefer in- 
teraction with a constant density 
medium 



e.g., 



Panaitescu Sz Kumar 



(I2OO2I ) who found that wind inter- 
action was preferred for just 1 burst 
(GRB 970508 ) out o f 10. However, 
Starling et al.l (120071 ) recently were 
able to constrain the circumburst 
medium for 5 BeppoSAX sources, 
finding that 4 were consistent with a 
wind and 1 (GRB 970508) was con- 
sistent with ISM. One difference with 



the an alysis of iPanaitescu fc Kumar 
(I2OO2I ) is that radio data were not 
included. In the Swift era, there have 
been excellent data on early X-ray 
afterglows, but there have been few 
extensive multiwavelength data sets. 
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A reason for this is that the greater 
sensitivity of Swift compared to the 
previous GRB satelhtes, so the bursts 
are fainter in multiwavelength ob- 
servations. In particular, there have 
been few radio hght curves, although 
the radio emission can provide im- 
portant constraints, as described 
above. 

Another aspect of Swift bursts is that 
the jet break often appears fairly late 
in the evolution, if at all. To some 
extent, this can be attributed to the 
discovery of lower luminosity bursts. 
An advantage of such bursts is that 
there is the possibility of using the 
distinguishing properties of wind vs. 
ISM models discussed above. A well- 
observed burst is GRB 050820A, 
which did not show a jet break un 



til an age >17 days (jCenko et al. 



20061 ) . During the pre-je t break pe- 
riod, ICenko et al.l (120061 ) found that 
the X-ray afterglow declines more 
rapidly than the optical afterglow, 
which is an indicator of ISM inter- 
action. However, an ISM model that 
is consistent with the optical and X- 
ray prope rties overpre d icts t he radio 



emission. iCenko et al.l (120061 ) expect 
a radio flux of 5 mJy on day 7, but 
observe a flux of 0.1 mJy. The low 
radio flux is consistent with a wind 
interaction model. Thus the situa- 
tion is ambiguous. 

In addition to the finding of late X- 
ray breaks in Swift bursts, optical ob- 
servations sometimes show a break 
when none i s present at X-ray wave - 
lengths (e.g. lMonfardini et al.l . l2006l ). 
This causes some uncertainty about 
the nature of the optical break and 
may indicate that the X-ray and op- 



tical emission come from different re- 
gions. 

Although multiwavelength model- 
ing provides the best constraints on 
afterglow models, the large set of X- 
ray light curves and spectra observed 
with Swift can be used for compari- 
son with the expected "closure rela- 
tions" for ISM and wind models. In 
a study of Swift bursts f rom the first 



6 mon ths of operation, IZhang et al. 



([20061) found that all the bursts were 
consistent with ISM interaction. 
However, when is below X-ray 
frequencies, the afterglow evolution 
does not depend on the density pro- 
file, limiting the number of objects 
for which an interesting result can be 
obtained. I n a more re c ent st udy of 
30 sources, iPanaitescul (120071 ) found 
that 2/3 are consistent with Uc below 
X-ray frequencies so the afterglow 
evolution does not depend on the 
density profile, 25% are consistent 
with ISM evolution, and 10% are 
consistent with wind evolution. 

The afterglow modeling described 
above applies to the blast wave phase 
of evolution in which the ejecta have 
been decelerated by the surrounding 
medium. The development of rapid 
response optical/infrared telescopes 
has given the possibility of making 
observations before the blast wave 
phase has been established. The 
REM telescope may have observed 
GRB 060418 and GRB 060607A dur- 
i ng the onset of the a fterglow phase 
(IMolinari et al.l . |2006| ) . Before decel- 
eration, the observed flux is expected 
to increase as t^ (' ISM) or t-^/^ ('wind) 



(iMohnari et all. l2006l : Ijin &: Fan 



20071 ) . The observed increases for the 
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2 bursts are consistent with a de- 
pendence and inconsistent with t^/^, 
implying a constant density interac- 
tion. An estimate of the radius at 
which decelerat i on oc curs is 10^^ cm 
( iMolinari et al.l . |2006| ) , showing that 
the constant density medium must 
extend in to at least this radius. An- 
other case where a sharp turn-on 
of the afterglow may have been de- 



tected is GRB 060206 f Stanek et al. 



20071 ). If this interpretation of the 
rise phases is correct, the further 
light curves should be of the ISM 
type (or something more complex), 
and not of the wind type. The avail- 
able information on these afterglows 
does not seem to fit the simple mod- 
els. Another burst with early optical 
observations is GR B 050801, which 
Rykoff et al.l (120061 ) found to have a 
flat flux evolution from 20 — 250 s. 
The early fiat evolution is roughly 
consistent with wind evolution, but 
the later afterglow evolution is con- 
sistent with ISM interaction, and not 
with wind interaction. It is possible 
that the burst made a transition from 
wind to constant density surrounding 
medium, but this would have to oc- 
cur close to the deceleration radius. 
Although there are tantalizing clues 
from the early optical observations, 
they cannot be clearly interpreted in 
terms of the standard models. 

One of the main findings during the 
Swift era is that X-ray afterglows are 
more complex than previously rec- 
ognized, showing a variety of flaring 
behavior and an early plateau phase. 
Detailed observations of optical af- 
terglows have al so shown c omple x 



evolution (e.g., iDai et al.l . 120071 ). 



These observations point to later en- 



ergy addition to the GRB blast wave 
than is assumed in the standard 
models. Our lack of knowledge of the 
expected form of the energy addition 
(unlike the supernova case) limits 
our ability to precisely deduce the 
nature of the surrounding medium. 
However, there is currently evidence 
for afterglow evolution in both wind 
and constant density media. 



6 Short vs. Long GRBs 



According to present thinking, the 
long GRBs are explosions in massive 
stars, while the short bursts result 
from the mergers of compact objects. 
These 2 progenitor types can be 
expected to have different environ- 
ments: the long bursts occurring in 
the mass loss of the progenitor stars 
and the short bursts in the surround- 
ing ISM. A comparison of the after- 
glows for the two types of bursts can 
then give an indication of whether 
the interpretation of long bursts in- 
teracting with a constant density 
medium is correct and is not the re- 
sult of an effect such as the variation 
of microphysical parameters. 

A good case is the analysis of GRB 
051221A bv ISoderberg et al.l (120061 ). 
The time of an apparent jet break 
was 5 days, so there was significant 
evolution in the pre-jet break regime. 
The X-ray afterglow decline was 
characterized by a = —1.06 ± 0.04, 
which, together with the X-ray spec- 
tral index, was consistent with evo- 
lution in the cooling regime. The 
flatter optical/X-ray spectral index 
and the flatter evolution at optical 
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wavelengths were roughly consistent 
with evolution in a constant density 
medium with z/c between optical and 
X-ray wavelengths, assuming a stan- 
dard afterglow model. This result 
gives confidence in the application 
of the standard model to those long 
bursts with the deduction that they 
are expanding into a constant density 
medium. The similarity between the 
short and long burst afterglow evolu- 
tion indicates that the apparent ISM 
interaction is not due to wind inter- 
action with a particular evolution of 
the microphysical parameters. 



The observations of the short bursts 
are generally consistent with interac- 
tion with a low densi ty ISM. In the 
caseo fGRB 05122lA. [So"derberg et al 
fl2006h deduced a density ~ 10"^ 



cm~^. The lack of observable X-ray 
afterglows for some short bursts may 
be due to a very lo w surroundin g 



density, <10-^ cm'^ flNakad . 120071 ). 



7 Producing a Constant Den- 
sity Surrounding Medium 



The evidence from afterglow mod- 
eling for constant density media 
around long GRBs has stimulated 
interest in producing such a medium 
around a massive star. The most 
plausible way of doing so is the 
medium produced downstream of 
the termination shock in the stellar 



wind (IWijersl . l200ll ) . This region has 
a roughly constant pressure because 
the sound speed is higher than the 
systematic velocities over most of 
the volume. In addition, the veloci- 
ties are sufficiently high to make the 



fiow steady over much of the volume, 
so that conservation of entropy with 
radius leads to a constant density re- 
gion. The radius of the termination 
shock, Rt, can be estimated from the 
pressure generated at the shock 



Rt = 5.7 X 10 
/ p/k 



19 



104 cm-3 K 



,103 km s-i 

-1/2 

AV' cm 



where p is the pressure in the shocked 
wind and k is Boltzmann's constant. 

A general problem is that the value 
of Rt needed to explain the afterglow 
observations is smaller than expected 
around a typical WR star. Afterglow 
models require that the termination 
shock be at a radius <2 x 10^'' cm in 



some cases ( IChevalier et al.l . l2004j ). 



One factor is the reduced value of 
M because of th e low rnetallic ity of 
the progenitor ( Wijerd . 2001 ). As 
discussed in section [2l this reduces 
M by a factor ~ 3, which reduces 
Rt by up to ~ 2. Although this 
helps the problem, more is needed. 
The other possibility is increasing 
the pressure, which can be accom- 
plished by interaction with a dense 
ambient medium, high ram pres- 
sure due to motion of the progenitor 
star, or high pressure of the ambi- 



ent medium (van Marie et al.l. 12006 



Chevaher et all. 2004). The finding 
of iFruchter et all fl2006l that GRBs 
occur in regions of strong star for- 
mation may be consistent with the 
presence of a high interstellar pres- 
sure, but this needs to be verified in 
more detail. 
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In the pre-burst models of lvan Marie et 
( I2OO6I ). the wind from the WR star 
sweeps out the dense red supergiant 
wind from a previous evolutionary 
phase. However, there is increasing 
evidence for supernovae occurring 
soon after the loss of the the H enve- 
lope in dense mass loss; an example 
is SN 2001em fIChugai fc Chevalieil . 
20061 ). GRBs might preferentially 
occur in such an explosion because 
there is less opportunity for loss of 
angular momentum in the WR wind. 
The dense mass loss then provides a 
wall for shocking the WR wind. 



A possible problem for the shocked 
wind explanation of the constant 
density medium is the radial range 
over which it is required. Observa- 
tions of some GRBs requires that 
the outer extent of the shocked wind 
be ^2Rt, which rules out some wind 



al.lThe evidence for the turn-on of some 
afterglows in an ISM medium (Sec- 
tion E]) is a problem for this model. 
The deceleration of the GRB ejecta 
in a wind occurs at a radius 



interaction models (j Chevalier et al. 
2OO4J ). While some shocked wind 
models are consistent with these re- 
sults, the lack of evidence for interac- 
tion with the region inside or outside 
the shocked WR wind is a possible 
problem. In particular, one would 
expect an interaction with a freely 
expanding wind followed by a tran- 
sition to constant density medium. 
Early work on this transition indi- 
cated that there would be increase 
in emission when the shock was tra- 



versed (IWiiersl.l200ll:lPe'er fc Wiier^ 



20061), but iNakar fc GranotI fl2006f ) 



find that there is no bump in the 
GRB light curve at this point. In any 
case, there should be a transition 
from the self-similar blast wave evo- 
lution in a wind medium to evolution 
in a constant density medium. 



R. 



dec 



4.0 X lO^^EsaFo 



cm, 



where £"53 is the isotropic blast 
in units of 10^^ 



wave energy 



ergs 



and Fq 2 is the initial Lorentz fac- 



tor o f the GRB ejecta in u nits o f 



10^ flPanaitescu fc Kumaii . [20001). 
The pressure needed to have the ter- 
mination shock occur at or within 
this ra dius cannot plausibly be at- 



tained (Ivan Marie et al.l . |2006| ). so a 



massive star progenitor may not be 
viable in these cases. A possibility 
is that some long bursts have com- 
pact binary progenitors a nd interact 



directly with the ISM. iKing et al 



( I2OO7I ) have suggested the merger 
of neutron stars and white dwarfs 
as possible long burst progenitors. 
These events would generally be as- 
sociated with active star formation 
(but not always) and would not be 
accompanied by a supernova. 



8 Clues from Absorption Lines 



Possible information on the imme- 
diate surroundings of GRBs comes 
from absorption lines observed in 
the optical spectra. The mass loss 
processes leading up to explosion 
have the possibility of creating ab- 
sorption features in the spectrum; 
in particular, the free wind is ex- 
pected to have a velocity as high 
as 5000 km s~^, which can be dis- 
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tinguished from typical velocities in 
the host galaxy. A redshift 2;>2 is 
needed so that the strong ultravi- 
olet resonance line transitions are 
shifted into the optical. A case of spe- 
cial interest has been GRB 021004 
[z = 2.3), which showed a num- 
ber of absorption line systems with 
velocities up t o 3000 km s~^ in ob - 
serv ations bv ISchaefer et al.l (120031 ) 



and iMirabal et aP fl2003h . who con 



eluded that the features were likely 
to be circumstellar. Higher resolu- 
tion observations have shown that 
the —3000 km s^^ line system can 
be separated into —2700 km s~^ an d 



—2900 km s ^ sy s tems ( iFiore et al. 



20051 : IChen et all . 120071 ). The prob- 



lem with the circumstellar hypothe- 
sis is that the radiation field of the 
GRB is expected to completely ion- 



stantial radius ( 


Lazzati et al.. 


2006; 


Chen et al.. 


2007 


)• 


Mirabal et al. 



( 120031 ) suggested that the high ve- 
locity features are due to clumps 
initially at some distance from the 
burst that are radiatively accelerated 
by the burst. However, it is not clear 
whether coherent acceleration by 
the radiation field could take place. 
An alternative point of view is that 
the 3000 km s~^ features are due to 
the freely expanding WR star wind, 
which has the advantage of naturally 
explaining the observed velocity. 
However, there is still t he problem of 



the st rong ionization. iLazzati et al 



( 120061 ) deal with this problem by 
suggesting that the WR wind termi- 
nation shock is out at 100 pc, where 
the free wind is not completely ion- 
ized. This value of Rt is in the op- 
posite direction to what is generally 
needed for GRB afterglows (Section 



[7]) and requires an unusually low sur- 
rounding pressure, in addition to the 
low surroun ding density a s sumed in 
the model of iLazzati et al.l (120061 ). 



Another problem with the WR star 
wind hypothesis is that the lines ob- 
served in spectrum of GRB 021004 
include H lines, which are generally 
not expected in WR star wind at the 
end of the star's life. The presence 
of H would not be a problem if the 
absorption i s formed in a n inte rven- 
ing system. IChen et al.l (120071 ) un- 
dertook a project to check whether 
absorption line systems in the ve- 
locity range 1000 — 5000 km s~^ 
were due to intervening systems or 
to the progenitor winds. The finding 
of one high velocity system out of 5 
observed GRBs was consistent with 
an intervening system. They also ar- 
gued that the high velocity systems 
observed in GRB 021004 could be 
attributed to an intervening system, 
citing the presence of H I, C II, and 
Si II together with the absence of ex- 
cited C II or Si II as evidence. Over- 
all, it appears that absorption line 
observations are not likely to give 
information on the medium that the 
GRB is exploding into. 

One way to demonstrate a rela- 
tion between the absorbing gas and 
the GRB would be line variabil- 
ity due to the GRB radiation. The 
C IV absorption features in GRB 
021004 were measured over 6 days 
but do not sho w any clear evidenc e 
for variability (Lazzati et al.l . |2006|). 
Vreeswiik etaO (l2006h found evi- 
dence for variability of Fe II and 
Ni II absorption lines toward GRB 
060418. They showed that ultraviolet 
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pumping of Fe II and Ni II excited- 
and metastable-level populations by 
the GRB radiation is plausible, but 
their model requires that the absorb- 
ing gas be ^1.7 kpc from the GRB. 
There is no significant Fe II or Ni II 
closer than this distance, presumably 
because of ionization by the GRB 
radiation. 



9 Discussion and Conclusions 



At the present time, there is a di- 
chotomy between the supernovae 
and the GRBs. In the supernova 
case, the driving force for the inter- 
action with surroundings is fairly 
well understood. The shock acceler- 
ation through the outer parts of a 
star does not depend on the details 
of the central explosion. Models for 
the nonrelativistic interaction with 
the surroundings of the progenitor 
star are straightforward and give 
strong support for models of interac- 
tion with a free wind. For GRBs, the 
relativistic, coUimated flow depends 
on the details of matter and energy 
production by the central engine. 
Observations of afterglows in the 
Swift era have shown that the ejecta 
properties are probably crucial for 
understanding the early afterglow 
evolution. 

The standard afterglow model is 
probably a better approximation at 
later times. From the beginning of 
afterglow observations, the standard 
interpretation of long burst after- 
glows have indicated that, in the 
majority of cases, the evolution is 
consistent with interaction with a 



constant densi t y med ium (except for 
Starling et all (l2007h l This trend 
has continued in the Swift era. The 
indications of a constant density 
medium include the evolution of the 
X-ray afterglow when the cooling 
frequency appears to be above X- 
ray wavelengths, and the evolution 
of the optical afterglow and its re- 
lation to the X-ray evolution. Radio 
observations do not support con- 
stant density interaction, and that 
is one of the main points against 
this picture. More recent support for 
constant density interaction in some 
cases comes from the sharp turn-on 
of optical afterglow emission. Also, 
the consistency of the standard con- 
stant density afterglow model with 
observations of short bursts gives 
confidence in the similar model when 
applied to long bursts; short bursts 
are expected to be interacting with a 
constant density ISM. In other cases, 
there is evidence for interaction with 
a wind medium. 

There is thus sufficient evidence for 
constant density interaction that 
ways of producing such a medium 
must be considered. The inner 
boundary of the constant density 
medium must extend in to <2 x 10^'' 
cm. This is larger than the radial dis- 
tance usually sampled by supernova 
observations, so there is not a clear 
discrepency between these cases. 
However, a constant density around 
a massive star is most plausibly pro- 
duced by having the stellar wind 
pass through a termination shock, 
and the required radius is smaller 
than would typically be expected for 
a massive star wind at the end of its 
life. There is thus interest in how to 
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produce a small termination shock 
radius for GRB progenitors. 

One effect is the lower mass loss 
rate expected for the low metallicity 
progenitors of the GRBs, but addi- 
tional effects are needed to bring in 
the radius. One possibility is that 
the GRBs occur in regions of the 
ISM with a high pressure, which re- 
duces the value of Rt. Some support 
for this is that the bursts are ob- 
served to be more concentrated to 
star fo rming regions than a re super- 
novae (IFruchter et al.l . l2006l ). The ac- 
tion of stellar winds and supernovae 
from massive stars can produce a 
high pressure. Other possibilities 
are a high space velocity in a dense 
medium or that the progenitor stars 
undergo a particular evolution that 
results in a small Rt- However, the 
observational evidence for interac- 
tion with a constant density medium 
at the GRB deceleration radius is 
difficult to reconcile with a shocked 
wind model and may indicate a com- 
pact binary progenitor in the ISM 
for these cases. In this situation, the 
GRB should not be accompanied by 
a supernova. The surroundings of 
GRBs may be indicative of various 
progenitor types. 
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